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' Abstract 

> ! 

■ Based on the extended optical model approach in which the polarization potential 

is decomposed into direct reaction (DR) and fusion parts, simultaneous \ 2 analyses are 
performed for elastic scattering, DR, and fusion cross section data for the 9 Be+ 208 Pb 
system at near-Coulomb-barrier energies. Similar \ 2 analyses are also performed by only 
taking into account the elastic scattering and fusion data as was previously done by the 
present authors, and the results are compared with those of the full analysis including 
the DR cross section data as well. We find that the analyses using only elastic scattering 
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and fusion data can produce very consistent and reliable predictions of cross sections 
particularly when the DR cross section data are not complete. Discussions are also given 
on the results obtained from similar analyses made earlier for the 9 Be+ 209 Bi system. 



PACS numbers : 24.10.-i, 25.70. Jj 
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I. INTRODUCTION 



D. 



In our recent study [1] , we have carried out simultaneous x 2 analyses of e 

"" ; 20SPK 2, 1 1 and 9 Be+ 209 Bi 



astic scatter- 



ing and fusion cross section data for the 6 Li+ 208 Pb XA y, U| and "tSe-K^tSi y, M systems 

□ n □ 

at near-Coulomb-barrier energies in the framework of an extended optical model [3, la |9( 
by introducing two types of complex polarization potentials: the direct reaction (DR) 
and fusion potentials. In such analyses, it is indispensable and essential to include the 
experimental data for the total DR cross section o e £ p and the fusion cross section a e p P , in 
addition to the elastic scattering cross section da e ^ p / dQ for the separate determination of 
the DR and fusion potentials. However, when the previous study was made Q, reliable 
data of a e ^ p for 6 Li+ 208 Pb and 9 Be+ 209 Bi were not available, and thus the analyses were 
proceeded in two steps. In the first step, we carried out \ 2 analyses of only the elastic 
scattering data by assuming just one simple Woods-Saxon type complex potential. Us- 
ing the potential parameters thus fixed, we could then generate the total reaction cross 
section or, which we called the semi-experimental total reaction cross section q-^" 11 -^ _ 
As has been shown in a number of publications, such ,j^ erm_e:r P predicted from the optical 
potential that fits the elastic scattering data usually reproduces a^ p very well. This is 

1(1 ll| , the deuteron 



the case for reactions induced by the proton 



and also heavy-ions 



121 , the a-particle 



LIS IS 



141 ]. We then used a R 



semi— exp 



to further extract semi-experimental 



i. i. l t~\t~> j_ • semi— exp 1 • ,1 i . • semi— exp semi— exp exp t 

total DR cross sections a D F by using the relation a D a R °f • m 

the second step, use was made of thus extracted a s ^ mi ~ ex P j n place of the experimental 

t^t~> j • exp ; j ■ i] l r i exp / in semi— exp j exp 

DR cross section a D to carry out simultaneous analyses of aa E /a\l } a D , and a F 
for determining the full extended optical model potential composed of two polarization 
potentials. 

The DR and fusion potentials thus determined revealed interesting characteristic fea- 



tures of these potentials. First of all, both potentials satisfy separately the dispersion re- 
lation Q| . Secondly, the fusion potential is found to exhibit a threshold anomaly 
as was observed for tightly bound projectiles 
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n 



15 



16] 



18l [l9| , but the DR potential does not 



is ueie 

as. 



show a pronounced threshold anomaly. Thirdly at the strong absorption radius, the 
magnitudes of the fusion potential were found to be much smaller than those of the DR 
potential. As a consequence, the resulting total polarization potential dominated by the 
DR potential becomes rather smooth as a function of the incident energy. This has solved 
a long standing puzzle why the threshold anomaly has not been seen in the polarization 
potentials_determined for the systems involving a loosely bound projectile such as 6 Li 
and 9 Be 

The extracted DR potentials have provided us with a unique opportunity to study 
the effects of breakup (DR) on fusion by comparing op calculated from either including 
or neglecting the real and the imaginary part of the DR potential. Such studies were 
made in Ref. 1J, which showed that in the sub-barrier region, breakup is not the main 
reason for the subbarrier enhancement of o e p V and that the mechanism that governs 
the enhancement is neutron flow as originally suggested by Stelson et al. 20]. In our 
approach, this effect is phenomenologically implemented in the imaginary part of the DR 
potential. On the other hand, in the above barrier region, the breakup suppresses a e p P 
and the observed suppression factors for 6 Li and 9 Be were fairly well accounted for in 
terms of the breakup 1|. 



After completing our work of Ref. ^| for 9 Be I -"'Bi system, clastic scattering data for 



3- 



9 Be + 208 Pb, a system similar to 9 Be + 209 Bi, have become available 211. Thus, for 9 Be 
_l_ sosp^ S ygt em we now have data available for the elastic scattering 2l| and fusion 
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cross sections as well as the sum of cross sections of breakup, transfer, and incomplete 
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fusion 23|. For a loosely bound projectile like 9 Be, we may assume the summed cross 
section thus observed as the total DR cross sections. This has provided us with an 
opportunity to carry out x 2 analyses taking into account all three sets of experimental 
data, i.e., the elastic scattering, DR, and fusion data. We can then compare these \ 2 
analysis results with those obtained by considering only two sets of data, i.e., the elastic 
scattering and fusion data without the experimental DR data. We shall henceforth call 
the case where all the three data sets are included in the x 2 analyses the EDF (elastic, 
DR, and fusion) approach, and the case where only two data sets are considered the EF 
(elastic scattering and fusion) approach. Note that in the EF approach we do, however, 
include a s £ m% ~ ex P^ wri ich is essential to fix the DR potential parameters. The aim of 
the present study is to make a comparison between these two approaches and study the 
validity of the EF method used in our previous work. By extending our EF method 
proposed in the previous work on 9 Be + 209 Bi to 9 Be + 208 Pb system, where we have DR 
data as well as elastic and fusion data, we shall show that the EF approach gives us very 
reliable predictions of cross sections. 

In Sec. II, we first generate a s ^ mi ~ ex P f or the EF approach case by following the method 
described in Ref. l|. Two types of x 2 analyses (EDF and EF) are then carried out in 
Sec. Ill and the results are compared and discussed in Sec. IV. Sec. V concludes the 
paper. 

II. EXTRACTING SEMI-EXPERIMENTAL DR CROSS SECTION 



Our method of generating a s ^ mi ~ ex P resorts to the empirical fact 24| that the total 
reaction cross section calculated from the optical model fit to the available elastic scat- 
tering cross section data, da e ^ p /dQ, usually agrees well with the experimental cxr, in 



spite of the well known ambiguities in the optical potential. Let us call the total reac- 
tion cross section thus generated the semi-experimental reaction cross section a s ^ mt - ex P_ 

rnu semi—exp „ . j i 

lhen, a D y is generated by 



semi—exp semi—exp exp 

a D -<t r -<t f . 



This approach seems to work even for loosely bound projectiles, as demonstrated recently 
by Kolata et al. [l4| for the 6 He+ 209 Bi system. 

Following Ref. l|, we first carry out rather simple optical model \ 2 analyses of elas- 
tic scattering data solely for the purpose of deducing a s ^ mi ~ ex P p or these preliminary 
analyses, we assume the optical potential to be a simple sum of two volume-type poten- 
tials Vo(r) and U\(r,E), where Vo(r) is the real, energy independent bare potential and 
Ui(r,E) is a complex potential with common geometrical parameters for both real and 
imaginary parts. The elastic scattering data are then fitted with a fixed radius parameter 
7*1 for U\(r,E) and with three other parameters adjustable; the real and the imaginary 
strengths V\ and W\ and the diffuseness parameter a±. The x 2 fitting is done for three 
choices of the radius parameter; ri—1.3, 1.4, and 1.5 fm. These different choices of the 
ri-value are made in order to examine the dependence of the resulting a s ^ m ' l ~ ex P on the 
value of 7*i. 

n 

As observed in Ref. |lj, the values of a s R emt exp thus extracted for three different re- 
values agree with the average within 1%, implying that ,j^ erm ~ e:r P is determined without 
much ambiguity. We then identified the average as the final value of aR emi ~ ex P Using thus 
determined a s R emi ~ ex P ) we generated a s ^ mi ~ ex P D y employing Eq. (1). The resultant values 

c semi—exp j semi—exp , j • rp ui t j .t. -.t. exp 

of a R and a D are presented m lable 1, together with a F 



23, a e * p 



j exp a c rr\ it t t\ i r semi—exp j semi—exp , . • 

and o R . As seen from lable 1, the values of a D and a R are systemati- 

cally larger than the corresponding experimental values, except for the lowest energy 



Table I: Measured and extracted fusion^DR, and total reaction cross sections for the 9 Be+ 208 Pb 
system. a £ p P and a e ^ p are from Refs. 



sion, UK, ana 1 

3. [2^1 and 2^|, respectively. a e £ p is the sum of a e p P and a e £ p . 



a 



semi— exp 
R 



is extracted from the elastic scattering data 



2 it ] as explained in the text. a D 



semi—exp 



is then obtained by using Eq. (1). 
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of E cm =36A MeV. The reason why a s ^ mi exp is larger than a e ^ p except for the lowest 
energy may be ascribed to the fact that a^ p includes contributions from only breakup, 



transfer, and incomplete-fusion events 



but not from inelastic scattering and other 



simple quasi-elastic processes such as pickups. The difference between a^ p and a s ^ mt ~ ex P 
becomes larger with energy. It implies that there are more open, but not identified, DR 
channels as the incident energy increases. 

It is worth remarking at this point that there is a reason to question the accuracy 
of the extracted value of a s ^ mt ~ ex P a t E cm =36A MeV. The experimental value of the 
ratio Pe of the elastic scattering to the Rutherford cross section at the forward angles 
are systematically larger than unity jsjj at this energy. The average value of Pe at 
small angles is about 1.033. This suggests that there may be a problem in the overall 



normalization constant in the measured data. In fact, it is indicated [2jj that there are 
experimental uncertainties of a few percents in the absolute normalization. Even just a 
few percent uncertainty in the normalization is critical, particularly at low energies in 
extracting cr ^ m4_ex P_ j n order to confirm this, we have reanalyzed the elastic scattering 
data by reducing the cross section by a factor of 1.033 so that the values of Pg at forward 
angles become around unity. A new value of a s ^ mi ~ ex P thus extracted turns out to be 
122 mb, which in turn gives us a s ^ mi ~ exp =112 mb. This value is significantly greater 
than 71 mb given in Table I and is also larger than the experimental value of a s ^ mt ~ ex P — 
109 mb. It is thus very plausible that the true values of a s ^ m ' l ~ ex P anc i ^mi-exp ^ 
energy could be larger than a e £ p and cr^ p , respectively. However, in the present x 2 
analyses, use is made of the a s ^ mz ~ ex P values as listed in Table I. 

III. SIMULTANEOUS x 2 ANALYSES 

Simultaneous x 2 ~ analyses are then performed for two cases of data sets; (da e F p /dQ, 

exp exp\ j I j exp /to semi—exp exp\ i j exp /to sxp j exp r ,1 

a D i a F ) and {d a F y aiL, a D , <J F ) where da F jd\l, a D , and a F are from the 



literatures 



21 



22 



231 ] . As mentioned in the Introduction, the former (latter) case with 
a D P (^a s ^ mi ~ exp ) is called the EDF (EF) analysis. In calculating the x 2 value, we simply 
assume 1% errors for all the experimental data. The 1% error is roughly the average of 
errors in the measured elastic scattering cross sections, but it is much smaller than the 
errors in the DR (~5%) and fusion (~10%) cross sections. The choice of the 1% error for 
DR and fusion cross sections is thus equivalent to increasing the weight for the DR and 
fusion cross sections in evaluating the x 2_va lues by factors of 25 and 100, respectively. 
Such a choice of errors may be reasonable, since we have only one datum point for each 
of these cross sections, while there are more than 50 data points for the elastic scattering 
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cross sections. 



A. Necessary Formulae 



The optical potential U(r, E) we use in the present work has the following form; 
U(r; E) = V c {r) - [V (r) + U F (r; E) + U D (r; E)], 



(2) 



where Vcif) is the usual Coulomb potential with r^=1.25 fm and Vq(t) is the bare 
(Hartree-Fock) nuclear potential. U F {r\ E) and Ud{t;E) are, respectively, fusion and 
DR parts of the polarization potential 2j| originating from couplings to the respective 
reaction channels. Both U F (r; E) and Ujj{r\ E) are complex and their forms are assumed 



to be of volume-type and surface-derivative-type 
and Uo(r; E) are explicitly given by 



2fij ]. respectively. Vo(r), Up(r;E), 



V (r) = V f(X ) 



(3) 



U F {r; E) = (V F (E) + iW F {E))f(X F ), (4) 

and 

U D (r; E) = (V D (E) + iW D (E))4a D ^^-, (5) 

where /(Xj) = [l+exp(Xj)] _1 with X{ = (r— Rj)/ai (i — 0, D and F) is the usual Woods- 
Saxon function, while V F (E), Vd(E), W f (E), and Wd(E) are the energy- dependent 
strength parameters. We assume the geometrical parameters of the real and imaginary 
potentials are the same, and thus the strength parameters Vi(E) and Wi(E) (i = F or 
D) are related through a dispersion relation 



where P stands for the principal value and Vi(E s ) is the value of Vi(E) at a reference 
energy E = E s . Later, we will use Eq. (6) to generate the final real strength parameters 
V F (E) and V D (E), after W F (E) and W D (E) are fixed from \ 2 analyses. Note that 
the breakup cross section may include contributions from both Coulomb and nuclear 
interactions, which implies that the direct reaction potential includes effects coming from 
not only the nuclear interaction, but also the Coulomb interaction. 

Vo(r) in Eq. (3) may also have an energy-dependence coming from the nonlocality 
due to the knockon-exchange contribution. We ignore such effects as they are expected 



to be small for heavy- ion scattering [27[ , and employ the real potential parameters used 
in Ref. assuming that all the unusual features of the potential may be put into 
the polarization parts, particularly in the DR part. The parameters used for Vq{t) are 
Vo=18.36 MeV, r =1.22 fm, and a =0.57 fm 3- 

Note that this potential is shallow, 
which is often required in fitting elastic scattering data of such projectiles as 6 Li and 



3 Be 



In performing the optical model calculation, one can evaluate o> and <je> by using the 



following expression 



IS 



30] 



— <x (+)\W l (r)\ X ^> 



For D) 



(7) 



where x is the usual distorted wave function that satisfies the Schrodinger equation 
with the full optical model potential U(r,E) in Eq. (2). o> and on are thus calculated 
within the same framework as da F /dQ is calculated. Such a unified description enables 
us to treat different types of reactions on the same footing. 
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B. Threshold Energies of Subbarrier Fusion and DR 



As in Ref. we also utilize as an important ingredient the so-called threshold ener- 
gies E ^f and £?o,c of subbarrier fusion and DR, respectively, which are defined as zero 
intercepts of the linear representation of the quantities Si(E), defined by 



Si = VEcri « an(E - E 0>i ) (i = DorF), (8) 

where is a constant. Si with i = F, i.e., Sp is the quantity introduced originally by 
Stelson et al. j^j], who showed that in the subbarrier region Sf from the measured op 
can be represented very well by a linear function of E (linear systematics) as in Eq. (8). 
In Ref. [26], we extended the linear systematics to DR cross sections. In fact the DR 
data are also well represented by a linear function. 

In Fig. 1(a), we present the experimental Sp(E) and Sp>(E). From the zeros of Si(E), 
one can deduce £ e J=30.0 MeV and E^J =35.0 MeV. For both i = F and D, the 
observed Si are very well approximated by straight lines in the subbarrier region and 
thus E 0ti can be extracted without much ambiguity. Another determination of Eq X £ can 
be made by using the semi-experimental DR cross section instead of the experimental 
DR cross section, as shown in Fig. 1(b). The resultant value, which we shall denote by 

E serni-e X p . g found ^ ^ ^emi-exp = 3^ ^ ^ ^ 

Eqj may then be used as the energy where the imaginary potential Wi(E) becomes 



zero, i.e., Wi(E 0ti ) = 



26 



. This procedure will be used later in obtaining a mathe- 



matical expression for Wi(E). 



11 



C. x 2 Analyses 



All the x 2 analyses performed in the present work are carried out by using Vq(t) as 
given in Subsec. Ill A and by using the fixed geometrical parameters for the polarization 
potentials, r>=1.40 fm, a^=0.30 fm, rn=1.50 fm, and ao=0.70 fm, which are close to 
the values used in our previous study [l|. Small changes of these values from the ones 



in 



used in Ref. [ll are made in order to improve the x 2 -fitt nig. 

As in Ref. |l|, the \ 2 analyses are done in two steps; in the first step, all 4 strength 
parameters, Vd(E), Wd(E), Vp(E) and Wp(E) are varied. In this step, we have been 
able to fix the strength parameters of the DR potential, Vu{E) and Wn(E), fairly well 
in the sense that the extracted Vd(E) and Wd(E) turn out to be smooth as functions of 
E. This is particularly the case for the imaginary strength Wd(E). The values of Vd{E) 
and Wd{E) are presented in Figs. 2 and 3 by open circles for the EDF and EF cases, 
respectively. It is remarkable that the resultant Wd(E) can be fairly well represented by 
the following function of E(= E c _ m _) (in units of MeV) 



W D (E) 



' for E < E e *l =30.0 

0.037(£ - 30.0) for 30.0< E <39.0 (9) 
0.33 for 39.0< E 



in the EDF case and 



W D (E) 



for E < E s ™- exp =32.5 

0.052(E - 32.5) for 32.5< E <40.0 (10) 

0.39 for 40.0< E 



in the EF case. Note that the threshold energies where Wjj(E) becomes zero are set equal 
to Eq X £ and E* e ™~ ex P as determined in the previous subsection and are indicated by the 
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open half circles sitting on the axis of E c/m . in Figs. 2 and 3. The dotted lines in the 
lower panels of Figs. 2 and 3 represent Eqs. (9) and (10), respectively. The dotted curves 
in the upper panels of Figs. 2 and 3 denote Vd as predicted by the dispersion relation 
Eq. (6), with W D (E) given by Eqs. (9) and (10), respectively. As seen, the dotted curves 
reproduce the open circles fairly well, indicating that Vd(E) and Wd(E) extracted by 
the x 2 analyses satisfy the dispersion relation. 

In this first step of x 2 fitting, however, V F (E) and W F (E) are not well fixed in the sense 
that the extracted values fluctuate considerably as functions of E. This is understandable 
from the expectation that the elastic scattering data can probe most accurately the optical 
potential in the peripheral region, which is nothing but the region characterized by the 
DR potential with r D = 1.5 fm. The part of the nuclear potential responsible for fusion 
with r F = 1.4 fm is thus difficult to pin down in this first step. 

In order to obtain more reliable information on V F and W F , we have thus performed 
the second step of the x 2 analysis. This time, instead of doing a 4-parameter search we 
use Vd and Wd determined by the first step of \ 2 fitting. But, rather than using Vd 
and Wd exactly as determined by the x 2 fitting, we use Wd(E) given by Eqs. (9) and 
(10) and V D (E) given by the dispersion relation. We then have performed 2-parameter 
X 2 analyses, treating only V F (E) and W F (E) as adjustable parameters. The values thus 
determined are presented in Figs. 2 and 3 by solid circles. As seen, both V F (E) and 
W F (E) are determined to be fairly smooth functions of E. The extracted W F (E) may 
be represented by 



W F (E) 



for E < E e J =35.0 

< 0.879(^-35.0) for 35. 0< E <38.3 



(11) 



2.90 



for 38.3< E 
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in the EDF case and 



W F (E) 







for E < =35.0 



0.771(^-35.0) for 35. 0< E <38.5 
2.70 for 38.5< E 



(12) 



in the EF case, respectively. As is done for Wd(E), the threshold energy where Wp(E) 
becomes zero is set equal to E^p and is indicated in Figs. 2 and 3 by the solid half circle 
on the axis of E c m . As seen, the Wp(E) values determined by the second x 2 analyses 
are fairly well represented by the functions given by Eqs.(ll) and (12). Note that the 



energy variations of Wp(E) and Vp(E) are quite rapid compared to those of 



IF 



17 



18, 



E) and 



191. 



Vd(E), and are similar to those observed in tightly bound projectiles 

Using Wp(E) given by Eqs. (11) and (12), one can generate Vp(E) from the dispersion 
relation. The results are shown by the solid curves in the upper panels of Figs. 2 and 3, 
which well reproduce the solid circles extracted from the x 2 -fhting. This means that the 
fusion potential determined from the present analysis satisfies the dispersion relation. 



D. Final Calculated Cross Sections in Comparison with the Data 

Using W D (E) given by Eqs. (9) and (10) and W F (E) given by Eqs. (11) and (12) 
together with Vd(E) and Vp(E) generated by the dispersion relation Eq. (6), we have 
performed the final calculations of the elastic, DR and fusion cross sections. Thus, instead 
of using the potential parameters just as extracted by the x 2 -analyses we have used these 
dispersive potentials for the final calculations. The results are presented in Figs. 4 and 
5 in comparison with the experimental data. All the data are well reproduced by the 
calculations, though there are subtle differences between the fits obtained by two types 
of the analyses as will be discussed in detail in Subsec. IV B. 
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IV. DISCUSSIONS 



A. Fusion and DR Potentials 

The characteristic features of the polarization potentials determined in the present \ 2 
analyses are very much similar to those obtained in our previous analyses jlj]. The real 
and imaginary parts of both fusion and DR potentials satisfy well the dispersion relation, 
and the fusion potential displays the threshold anomaly. As already presented in Figs. 
2 and 3, these features are seen in the strength parameters, Vp(E), Wp(E), Vd(E) and 
W D (E). 

Another important feature of the extracted potentials is that at the strong absorption 
radius of R sa =12.3 fm both the real and imaginary parts of the DR potential are consider- 
ably greater than those of the fusion potential, although the strength parameters Vd(E) 
and Wd(E) are smaller than Vp(E) and Wp(E). Thus, the energy dependence of the net 
polarization potential (sum of the fusion and DR potentials) at R sa becomes dominated 
by the DR potential which has a relatively smooth energy dependency. Consequently, 



the net potential does not show such a thresho 



for systems with tightly bound projectiles 



17 



d anomaly as seen in the net potential 



18, 



3 



However, after separating the 



polarization potential into DR and fusion parts, we clearly observe the characteristic 
threshold anomaly in the fusion potential. 



B. Comparison of EDF and EF Cross Sections 

Both EDF and EF approaches reproduce the experimental elastic scattering cross 
sections very well as shown in Fig. 4. The calculated cross sections shown in the left and 
right panels of Fig. 4 agree well with each other. It may then be naturally expected that 
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the resultant total reaction cross sections also agree with each other. This is indeed the 
case; the values of the calculated total reaction cross sections from the EDF approach 
are approximately equal to those from the EF approach, as shown by the dashed curves 
in Figs. 5(a) and (b). 

Since a s ^ mi ~ ex P j s extracted from the fit to the elastic scattering data, our final cal- 
culation results using the dispersive potential naturally reproduce a s ^ mt ~ ex P as shown by 
the dotted curves in Fig. 5(b). In the EF case, the calculations also reproduce both DR 
and fusion cross sections as well. This is, however, not the case for the EDF approach; 
calculations using the dispersive potential somewhat overestimates the experimental data 
of all three cross sections as Fig. 5(a) shows. 

It may thus be concluded that the overall fit to the data obtained in the EF case is 
better than the EDF case and that the main source of problems in getting a goo d overall 
fit in the EDF case comes from inconsistency between the elastic scattering j2jj and the 
DR reaction 2^ data; the elastic scattering data require more absorption (larger total 
reaction cross section) than what the measured total absorption (reaction) cross sections 
tell us. In view of this, it is important that measurements be made of inelastic scattering 
and some other quasi-elastic reactions which are not taken into account in the total DR 
cross section used in the present analyses. 

C. Effects of Breakup on Fusion 

We now turn to the effect of breakup on the fusion cross section. As has been argued, 
there are two competing physical effects of breakup on the fusion cross section, op. The 
first is the lowering of the fusion barrier, which tends to enhance op. The other is the 
removal of flux from the elastic into the breakup channel, which suppresses <jp. Since the 

16 



breakup channel dominates DR, these two competing breakup effects may be represented 
by the real (Vb(r; E)) and the imaginary (Wp>(r; E)) parts of the DR potential; Vp(r, E) 
can describe precisely the effect of lowering the barrier, while Wp(r,E) the removal of 
the flux from the elastic channel. 

To see the effects quantitatively, we have introduced in Ref. [lj] the following suppres- 
sion factor R th , 

R* = a F /a F (V D = W D = 0), (13) 

where (Tp(Vp = Wp — 0) is o> obtained by setting Vp = Wp = 0, i.e., neglecting both 
barrier-lowering and flux loss effects, while op is our final calculated cross section that 
includes both Vp and Wp. In the above-barrier region, R th becomes almost constant and 
here we present just the average of the R th - values at three highest energies considered in 
the present study. The values are 0.87 and 0.82 for the EDF and EF cases, respectively. 
Setting Vp — reduces a F , while setting Wp = increases op. Thus, the fact that the 
.Revalues are smaller than unity indicates that the flux loss effect surpasses the barrier- 
lowering effect in the above barrier region. The theoretical values may be compared with 
the experimental values of R exp =0.79, where R exp is defined as 

R exp = <Tf"/<t f {V d = W D = 0). (14) 

by using <7p(Vp = Wp = 0) fixed from the EF case. R exp in the EDF case is 0.77, quite 
close to R exp in the EF case. 

Note that the theoretical suppression factor R th =0.82 in the case of EF agrees very 
well with the experimental value of R exp =0.79. It is natural because the calculated op 
agrees with a e p V in the EF shown in Fig. 5(b). Similarly, the difference between 

R th =0.87 and R exp =0.77 in the EDF case originates from the discrepancy between <j f xv 
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and the calculated op, seen in Fig. 5(a). In either case, both R th and R exp are consistently 
and considerably smaller than unity, implying that the observed suppression of op can 
be ascribed to the flux loss in the elastic channel to breakup. A similar result was also 
obtained in Ref. 1|. 

Although breakup (or DR) is the dominant factor in the suppression of op in the 
above barrier region, this is not the case in the sub-barrier region, where the neutron flow 
affects fusion dominantly 2(|, generally enhancing the sub-barrier fusion. In Ref. jl|, it 
was proposed that a good measure for the sub-barrier fusion enhancement is the quantity 
A defined as 

A = V B - Eojr, (15) 

where Vb is the Coulomb-barrier height and Eo t p is the sub-barrier threshold energy 
discussed in Subsec. II B. In Ref. lj, it is demonstrated that A is very well proportional 
to the neutron transfer Q- value. 



D. Comments on the Analyses of the 9 Be+ 209 Bi System Reported in Ref. [1 



In Ref. 1|, we presented our analyses on the 9 Be+ 209 Bi system using only the elastic 



scattering 



5J] and fusion cross section data p (the EF type analysis). Since the target 



nucleus 209 Bi differ from 208 Pb only by one proton, it is naturally expected that the 
experimental cross sections for the two systems should be very similar. This is indeed 
the case for the elastic scattering cross sections; no noticeable difference can be found 



in the data measured for the Pb target 



2l| and Bi target 



a. 



In contrast to this, the 



values of the fusion cross section for the Bi target we used from Ref. 



at the time of 



our analyses |l| are significantly larger than those for the Pb target reported in Ref. 
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Recently, however, the fusion cross sections for the Bi target are revised |32|, and the 
revised values are now very much the same as those of the Pb target. 

Due to this change in the experimental values of a e p P for 209 Bi, we have repeated our 
previous analyses for 9 Be + 209 Bi system, obtaining now essentially the same results as 
in the present work for 9 Be + 208 Pb. Therefore, we take this opportunity to revise our 
previous values of the suppression factor R; the new theoretical value obtained with the 
revised data is R th =0.81, which can be compared with the new experimental value of 
R exp =0.79. The corresponding values reported previously in Ref. jl| were R th =0.89 and 
R ex P=0.92. 



V. CONCLUSIONS 

In summary, we have carried out simultaneous \ 2 analyses of elastic scattering, DR 
(breakup plus incomplete fusion), and fusion cross sections for the 9 Be+ 208 Pb system 
at near- Coulomb-barrier energies within the framework of an extended optical model 
that introduces the DR and fusion potentials. Two types of analyses are made; one 
using the experimental DR cross section <jp P (EDF case), and the other using the semi- 
experimental DR cross section a s ^ mi ~ ex P (EF case), together with the measured elastic 
scattering and fusion cross sections for both cases. In the second type of the analyses, 
^semi-exp - g g rg ^. ex £ rac t ec [ f rom simple optical model fits to the elastic scattering data 
only. The extracted a s ^ mi ~ ex P are found to be significantly larger than a e ^ p . In spite of 
this difference between a e £ p and a s ^ mi ~ ex P ; the resultant DR and fusion potentials show 
common features that they satisfy fairly well the dispersion relation and the fusion 
potentials show the threshold anomaly as seen in the potentials for systems with tightly 

,3- 



bound projectiles 
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17. 



18, 
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For both EDF and EF cases the elastic scattering cross sections are equally well 
reproduced. However, the calculated DR, fusion, and total reaction cross sections fit the 
corresponding experimental data well in the EF case, but not in the EDF case. In the 
latter case, the calculations overestimate significantly the experimental DR, fusion, and 
total reaction cross sections. This is because there are some reaction channels that are 
not taken into consideration in the present experimental DR data. 

Thus as far as we don't have comprehensive a^ p available, the EF analysis gives us 
better overall results than the EDF analysis. We believe that if the cross section of 
inelastic scattering and some other missing reactions that are not taken into account 
in the present data of a e ^ p are measured and used in the analyses, both types of 
analyses will lead to equally good fit to the data. It is thus highly desirable that such 
DR data will be measured in near future in order to test our expectation and thus to 
justify the validity of the EF method proposed in Ref. 
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E c. m . < MeV ) 



Fig. 1: The Stelson plot of Si = y/E c , m ai for DR (i = D, open circles) and fusion (i = F, 
solid circles) cross sections with (a) the experimental and (b) the semi-experimental DR cross 
sections. The straight lines are drawn to show the extraction of the threshold energies Eq^. 
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Fig. 2: The strength parameters Vi (upper panel) and Wi (lower panel) for i = D and F as 
functions of £7 c . m . in the EDF case. The open and solid circles are the strength parameters 
for % = D and F, respectively. The dotted and solid lines in the lower panel denote Wd and 
Wf from Eqs. (9) and (11), respectively, while the dotted and solid curves in the upper panel 
represent Vd and Vf calculated by using the dispersion relation of Eq. (6) with W{ given by 
Eqs. (9) and (11). The reference energies, Vf(E s ) and Vd(E s ), are chosen as 4.0 MeV and 
0.85 MeV, respectively. 
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Fig. 3: The same as in Fig. 2, but for the EF case. The dotted (solid) line in the lower 
panel denotes Wr> (Wf) from Eq. (10) (Eq. (12)). The dotted (solid) curve in the upper panel 
represents Vd (Vf) obtained by the dispersion relation. The reference energy E a for Vf(E s ) 
and Vd(E s ), are taken as 3.2 MeV and 0.85 MeV, respectively. 
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Fig. 4: Ratios of the elastic scattering cross sections to the Rutherford cross section calculated 
with our final dispersive optical potential for (a) the EDF and (b) the EF cases are shown in 
comparison with the experimental data. The data are taken from Ref. J^. 
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Fig. 5: DR and fusion cross sections calculated with our final dispersive optical potential for 
the (a) EDF and (b) EF cases are shown in comparison with the experimental data. a s ^ irit ~ ex P 
denoted by the open circles in EF case are obtained as described in Sec. II. a e ^ p denoted by the 

n 

open circles in EDF case are the experimental DR cross sections [23|]. The fusion data are from 
Ref. Q 
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FIGURE CAPTIONS 



Fig. 1. The Stelson plot of Si = \/E c . m Oi for DR (i = D, open circles) and fusion (i = 
F, solid circles) cross sections with (a) the experimental and (b) the semi-experimental 
DR cross sections. The straight lines are drawn to show the extraction of the threshold 
energies E 0ti . 

Fig. 2. The strength parameters Vi (upper panel) and Wi (lower panel) for i = D and 
F as functions of E c-m , in the EDF case. The open and solid circles are the strength 
parameters extracted by x 2 -htting for % = D and F, respectively. The dotted and solid 
lines in the lower panel denote Wd and Wp from Eqs. (9) and (11), respectively, while 
the dotted and solid curves in the upper panel represent Vd and Vf calculated by using 
the dispersion relation of Eq. (6) with Wi given by Eqs. (9) and (11). The reference 
energy E s for Vp(E s ) and Vd(E s ), are chosen as 4.0 MeV and 0.85 MeV, respectively. 
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Fig. 3. The same as in Fig. 2, but for the EF case. The dotted (solid) line in the lower 
panel denotes Wd (Wf) from Eq. (10) (Eq. (12)). The dotted (solid) curve in the upper 
panel represents (Vp) obtained by the dispersion relation. The reference energy E s 
for Vf(E s ) and Vd(E s ), are taken as 3.2 MeV and 0.85 MeV, respectively. 



Fig. 4. Ratios of the elastic scattering cross sections to the Rutherford cross section 
calculated with our final dispersive optical potential for (a) the EDF and (b) the EF cases 
are shown in comparison with the experimental data. The data are taken from Ref. (21 ]. 



Fig. 5. DR and fusion cross sections calculated with our final dispersive optical poten- 
tial for the (a) EDF and (b) EF cases are shown in comparison with the experimental 
data. a s 'r m ' l ~ ex P denoted by the open circles in EF case are obtained as described in Sec. II. 

n 

a D denoted by the open circles in EDF case are the experimental DR cross sections [23]. 
The fusion data are from Ref. 1221. 
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